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ABSTRACT 

We present the third data release from the Australia Telescope Large Area Survey 
(ATLAS). These data combine the observations at 1.4 GHz before and after upgrades 
to the Australia Telescope Compact Array reaching a sensitivity of 14p,Jy beam -1 
in 3.6 deg 2 over the Chandra Deep Field South (CDFS) and of 17 /iJy beam -1 in 
2.7 deg 2 over the European Large Area ISO Survey South 1 (ELAIS-S1). We used a 
variety of array configurations to maximise the uv coverage resulting in a resolution of 
16 by 7 arcsec in CDFS and of 12 by 8 arcsec in ELAIS-S1. After correcting for peak 
bias and bandwidth smearing, we find a total of 3034 radio source components above 
5cr in CDFS, of which 514 (17 per cent) are considered to be extended. The number 
of components detected above 5er in ELAIS-S1 is 2084, of which 392 (19 per cent) 
are classified as extended. The catalogues include reliable spectral indices (Aa < 0.2) 
between 1.40 and 1.71 GHz for ~ 350 of the brightest components. 

Key words: catalogues — radio continuum: galaxies — surveys — methods: data 
analysis 


1 INTRODUCTION 

Large multiwavelength surveys are indispensable for under¬ 
standing galaxy formation and evolution. Radio wavelengths 
are valuable in providing an obscuration-independent tracer 
of star formation and active galactic nucleus (AGN) activity. 
However, most radio surveys so far have either covered only 
small areas, thus suffering from sample and cosmic variance, 
and missing intrinsically rare objects, or cover wide areas but 
are relatively insensitive, and therefore miss the most active 
epochs of galaxy formation. In particular, most wide radio 
surveys have not had sufficient sensitivity to detect normal 
star formation activity in any but the most nearby galax¬ 
ies, limiting their ability to contribute to our understanding 

*Email: thomas.franzen@curtin.edu.au 


of the cosmic evolution of galaxies. Another important req¬ 
uisite to maximising the astrophysical value of a survey at 
any wavelength is to maximise the overlap with other wave¬ 
lengths. 

Here we present the third data release of the Aus¬ 
tralia Telescope Large Area Survey (ATLAS), which cov¬ 
ers ~ 6 deg 2 to an rms depth of ~ 15/tjy beam -1 . The 
ATLAS survey area consists of t wo regions centred on the 
Chandra Deep Field South ('CDFS: fGiacconi et al.ll200llf and 
the European Lar ge Area ISO Survey - South 1 (ELAIS-S1; 
I Oliver et al. im. These two fields were carefully chosen 
to coincide with areas imaged by the Spitzer Wide-area In - 
frared Extragalactic Survey fSWIRE: iLonsdale et al]|2003h 
program, so that infrared and optical data are available for 
most of the radio objects. They also encompass the well- 
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studied Great Observatories Origins D eep Survey (GOODS) 
field in CDFS dGiavalisco et al.1 2004) 1. 

These two areas have since been the target of 
many other deep multi-wavelength sur veys, such as th e 
4 Ms X-ray survey from Chandra dXue et al.l 1201 ill . 
the Spit zer Extragalactic R epresentative Volume Survey 
(SERVS: lMauduit et al.ll2012ll. the VISTA Dee p Extragalac¬ 
tic Observations ( VIDEO: I Jarvis et aid l2013l l survey and 
the Herschel Multi-tiered Extragalactic Survey (HERMES; 
lOliver et al.ll2012h as well as spectroscopic and photometric 
redshift su rveys, including t he PRIsm MUlt i-object Survey 
(PRIMUS; ICoil et al.l[201ll : ICool et aldEoill. the FourSta r 
Galaxy Evolution Survey (zFOU RGE; ISpi tler et aI1l2012lh 
the ATLAS spectroscopy program dMao et alj|2012li ~and the 
new OzDES spectroscopy program (Lidman et al., in prepa¬ 
ration), making them some of the best-studied fields in the 
sky. 

The first ATLAS data release (DR1: [Norris et al]l2006l : 
iMiddelberg et akll2008l l surveyed these two fields in CDFS 
and ELAIS-S1 to a typical rms of 36 and 29 pj y beam - 1 
respec tively. The second data release (DR2: lHales et aid 
l2014al lbl) surveyed them to a typical rms of 30 /rJy beam -1 , 
addressing a number of sources of systematic error in DR1 
and for the first time presenting ATLAS polarisation results. 
Here in DR3 we report a further improvement to respective 
rms sensitivities of 14 and 17 pJy beam -1 . 

Minimising the rms noise is critical, because deep radio 
surveys such as ATLAS probe flux densities approaching the 
point where star-forming galaxies start to dominate the ra¬ 
dio sky. Surveys of radio sources with flux densities greater 
than 1 mjy are typically dominated by AGNs, but sources 
at lower flux dens ities are increasingly driven by star forma¬ 
tion activity (e.g. Seymour et al.l 120081 : ISmolcic et al.l 120081 : 
iBonzini et a,lj|2013 1.' Since star-forming galaxies also dom¬ 
inate non-radio surveys, the fraction of optical/IR galaxies 
detected at radio wavelengths rises sharply with decreasing 
flux density, yielding an extinction-free measure of star for¬ 
mation rate. 

The key science goals of ATLAS are: 


• To determine the relative contribution of starbursts and 
AGN to the overall energy density of the universe, and the 
relationship between AGN and star-forming activity. Par¬ 
ticularly interesting are composite galaxies in which a radio 
AGN lies buried within a host galaxy whose optical/infrared 
spectrum or SED appears to be that of a star-forming 
galaxy. 

• To test whether the far-infrared-radio correlation 
changes with redshift or with other galaxy properties. Once 
calibrated, this correlation will be an effective method to 
measure the star formation history of the Universe. 

• To trace the radio luminosity function to a high redshift 
(z ~ 1) for moderate-power sources, and measure the differ¬ 
ential 20 cm source count to a flux density limit of ~ 100 p Jy 
to a high precision. 

• To explore regions of parameter space, corresponding 
to a large area of sky surveyed at high sensitivity at radio, 
mid-infrared, and far-infrared wavelengths, to discover rare 
but important objects, such as short-lived phases in galaxy 
evolution. 


A supplementary but important goal is to act as 
a pathfinder for the Evolutionary Map of the Universe 


(EMU; iNorris et al.l 1201 il l survey, w hich will use the new 
Australian SKA Pathfinder (ASKAP: Ijohnston et al1l2007l . 
120081 : iDeBoer et al.ll2009h telescope to make a deep (10/rJy 
beam -1 rms) radio continuum survey of the entire Southern 
Sky, extending as far North as +30°. EMU will cover roughly 
the sa me fraction (75%) of the sky as the benchmark NVSS 
survey dCondon et al.lll998l l , but will be 45 times more sensi¬ 
tive, and will have an angular resolution (10 arcsec) 4.5 times 
higher. EMU is expected to generate a catalogue of about 70 
million galaxies, compared to the ~ 2.5 million sources cur¬ 
rently known at all radio frequencies. Since EMU will have 
a similar resolution to ATLAS, we are using ATLAS to test 
many of the technical and scientific processes for EMU. In 
addition, we will use optical spectroscopy of ATLAS galaxies 
to train the photometric redshift algorithms for EMU. 

Throughout this paper we define a radio ‘component’ 
as a discrete region of radio emission identified in the source 
extraction process. We define a radio ‘source’ as one or more 
radio components that appear to be physically connected to 
one host galaxy. Thus, we count a classical triple radio-loud 
source as being a radio source consisting of three radio com¬ 
ponents, but count a pair of interacting starburst galaxies 
as being two sources, each with one radio component (pro¬ 
viding of course that the angular separation between the 
starburst galaxies is large enough for them to be resolved). 

This paper (Paper I) is primarily concerned with de¬ 
scribing the survey and presenting the component catalogue. 
Paper II (Banfield et al. 2015, in preparation) will extract 
the component counts and explore the distribution of spec¬ 
tral indices, and Paper III (Norris et al. 2015, in preparation) 
will group the components into sources with optical/infrared 
identifications. Subsequent papers in the series will address 
the individual science goals of ATLAS. 

This paper is organised as follows. Section [2] describes 
the observations, calibration and imaging, and Section [3] de¬ 
tails the corrections made to the final image prior to com¬ 
ponent extraction. Section [4] describes the component ex¬ 
traction process and presents the ATLAS DR3 component 
catalogue. We summarise our results in Section [5] 

Throughout this paper we assume a Hubble constant 
of 71 km s -1 Mpc -1 , and matter and cosmological constant 
density parameters of Dm = 0.27 and Da = 0.73, and as¬ 
sume the convention for spectral index, a, where S oc v a . 


2 OBSERVATIONS, CALIBRATION AND 
IMAGING 

2.1 Target Fields 

ATLAS covers two regions, each of ~ 3 deg 2 , surrounding 
CDFS (a = 03 h 30 m 16.3 s , 8 = -28 o 05'12.4") and ELAIS-S1 
(a = 00 h 33 m 50.8 s , <5 = — 43°44'57.4"). These two fields have 
previously been observed at 1 .4 GHz as part of the ATLAS 
project and are described by lHales et al.l (l2014a| j. Both of 
these fields were originally targeted for radio observations as 
they overlap with SWIRE which includes infrared and opti¬ 
cal data for the majority of the radio objects. Tabic Q] lists 
the available radio observations (74 MHz - 20 GHz) overlap¬ 
ping the two ATLAS fields. 

To compare ATLAS with other existing 1.4 GHz surveys 
shown in Fig. [l] we compile a figure of merit for surveys, 
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Table 1. Available radio data and observations in CDFS and ELAIS-S1. Listed is the wavelength, telescope and survey where applicable, 
area overlap with the data presented in this paper, sensitivity, angular resolution and reference. 


Field 

Wavelength 

Survey/ 

Overlapping 

Sensitivity 

Resolution 

Reference 



Instrument 

Area (deg^) 

(MJy/bm) 

( arcsec) 


CDFS 

15 mm (20 GHz) 

ATCA 

2.50 

300 

29.1 X 21.9 

Franzen et a^^20^4) 


6 cm (5.5 GHz) 

ATLAS/ATCA 

0.25 

12 

4.0 X 2.0 

Huvnh et a^^^0^2) 


13 cm (2.3 GHz) 

ATLAS/ATCA 

3.57 

80 

57.2 x 22.7 

Zinn eiwi^^^O^) 


21cm (1.4 GHz) 

VLA 

0.33 

7.4 

2.8 X 1.6 

Miller et al. ('20131 


21 cm (1.4 GHz) 

ATLAS/ATCA 

3.57 

40 

11 X 5 

Norris et al.^00£) 


21 cm (1.4 GHz) 

ATLAS/ATCA 

3.57 

30 

12 X 6 

Hales et al. (2014a) 


21 cm (1.4 GHz) 

NVSS/VLA 

3.57 

450 

45 

Condon et a^^^9^1 


92 cm (325 MHz)* 

GMRT 

3.57 

100 

7.1 x 10.8 

Afonso et al. (^007) 


1.25 - 3.75 m (80 - 240 MHz)* 

GLEAM 

3.57 

3000 

150 

Wavth et al. (^0H5) 


1.5 - 3m (104 - 196 MHz) 

MWACS 

3.57 

40000 

180 

Hurlev-Walker et al. (2014) 


2 m (150 MHz) 

TGSS/GMRT 

3.57 

5 — 7 mjy 

20 

Sirothia et a^^^O^) 


4 m (74 MHz) 

VLSS/VLA 

3.57 

0.1 Jy 

80 

Cohen et al. (2007) 

ELAIS-Sl 

13 cm (2.3 GHz) 

ATLAS/ATCA 

2.70 

70 

33.6 x 19.9 

Zinn et al. ^0^) 


21 cm (1.4 GHz) 

ATLAS/ATCA 

2.70 

30 

10.3 x 7.2 

Middelberg et al. (2008) 


21 cm (1.4 GHz) 

ATLAS/ATCA 

2.70 

30 

12 x 6 

Hales et al. (2014a) 


36 cm (843 MHz) 

MOST 

2.70 

600 

62 x 43 

Randall et a^^20^2) 


50 cm (610 MHz) 

GMRT 

3.30 

100 

10 X 5 

Interna I^^M) 


2 m (150 MHz) 

TGSS/GMRT 

2.70 

5 — 7 mjy 

20 

Sirothia et a^^20^) 


1.25 - 3.75m (80 - 240 MHz)* 

GLEAM 

2.70 

3000 

150 

Wavth et al. (^ORT) 


1.5 - 3m (104 - 196 MHz) 

MWACS 

2.70 

40000 

180 

Hurlev-Walker et al. (2014) 


* Data yet to be released. 


SfoM ■ The goal of a survey is to maximise the area observed 
and minimise the noise in the image. This is limited by the 
fact that it takes T 2 times longer to decrease the thermal 
noise by a factor of T and T times longer to increase the 
observing area by a factor of T. Therefore, Sf 0 m is of the 
form: 

SFoM= (^’ (1) 

where S7 is the survey area in square degr ees and 5n m is 
the 5c r detection limit of the survey in mjy. iBunton fe Havl 
(l20ldl used a similar metric to quantify the survey speed for 
chequerboard phased array feeds. Of the current 1.4 GHz 
surveys, FIRST has the largest Sf 0 m (1 x 10 4 deg 2 mjy -2 ). 
ATLAS has a similar SfoM (1108 deg 2 mjy -2 ) to LH-Owen, 
COSMOS and Stripe82, thereby providing among the deep¬ 
est and widest coverage of radio objects at 1.4 GHz. 

We used the Australia Telescope C ompact Array 
(ATCA; iFrater, Brooks, fe Whiteoakl Il992l l at 1.4 GHz to 
observe the two ATLAS fields in all four Stokes parame¬ 
ters (XX,YY,XY,YX). During our observing campaign, 
the ATCA was being upgra ded with the Comp act Array 
Broadband Backend (CABB; IWilson et al.ll201lf ) providing 
a larger instantaneous bandwidth coupled with increased 
sensitivity of continuum and spectral line observations. At 
the time of our ATLAS observations, the CABB band pro¬ 
vided a 500 MHz bandwidth covering 1.3 — 1.8 GHz, split 
into 1 MHz channels. 

ATCA project C196T was allocated 1000 hours dis¬ 
tributed over 78 days between 2009 June and 2010 June to 
extend ATLAS. The observations were spread over the four 
6 km array configurations to maximise uv coverage as listed 
in Table [XT1 The primary flux a nd bandpass cal ibrator PKS 
1934-638 (14.95 Jy at 1.380 GHz: lRevnoldslll994l j was visited 
at the beginning of each observing run. Two different sec¬ 
ondary phase calibrators, PKS 0022-433 for ELAIS-S1 and 
PKS 0400-319 for CDFS, were observed every 15 minutes 
throughout the observations to determine the antenna com¬ 
plex gains and polarization leakage correction. We observed 


the two fields in the standard ATCA mosaic mode with 28 
pointings in CDFS and 20 pointings in ELAIS-S1. The 48 
pointing centres are listed in Tabic f/Vll 

2.2 Calibration 

We calibrated and edited th e uv data using MIRIAlQ 
dSault. Teuben, fe Wrightlll995h . The standard MIRIAD cali¬ 
bration techniques are optimised for the original ATCA cor¬ 
relator system (2 x 128 MHz bandwidths) and had to be 
expanded to calibrate the new CABB-enabled bandwidth. 
We calibrated each day of observations separately following 
the method outlined below. 

We removed channels affected by self-interference due 
to the 640-MHz clock harmonics and by known radio fre¬ 
quency interference with atlod in MIRIAD using the options 
birdie and rfiflag. We then restricted the frequencies to 
the known range of good bandpass response using uvaver 
to include only data between 1.3 and 1.8 GHz. Calibration 
of the bandpass was completed using the standard MIRIAD 
procedure mfcal and we applied the calibration to the sec¬ 
ondary calibrator. Automatic flagging of the calibrators was 
completed using mirf lag, setting int=600. This task imple- 
ments the automa tic flagging routine pief lag developed by 
iMiddelberd d2006h . We also manually flagged the calibrators 
using uvflag. 

The large bandwidth poses an issue with frequency de¬ 
pendent calibration. The standard MIRIAD calibration pro¬ 
cedure was altered to deal with the large bandwidth when 
calibrating the antenna gains and phases and the instru¬ 
mental polarization. The task gpcal has been updated to 
allow for calibration over smaller frequency bands instead of 
over the full CABB band. For our data, we calibrated over 
128 MHz sections across the band using the option nbins 
in gpcal. Once the calibration was complete, we copied the 
solutions to the targets using gpcopy. The target pointings 

1 http://www. atnf.csiro.au/computing/software/miriad/ 
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Figure 1. Comparison of ATLAS with the following 1.4 GHz sur¬ 
veys: the Evolutionary Map of the Universe fEMU ^Norris et al.l 
1 20 111) survey, the Westerbork O bservations of th e Deep APER- 
TIF Northern-Sky (WODAN; iRottgeri ng et al ] 120101) survey, 
the NR AO VLA Sky Survey (NVSS: ICondon et al.l 1 1998h . the 


iBecker. White fe He 

fane 

19951) survey, the Phoenix Deep Sur- 

vey I Hopkins et ak[ 

2003 

), the Australia Telescope ESO Slice 


the Bootes field by_ de Vries et al.l(| 20021 ) , th e Marano Field by 
iGruppioni et all (|l997l ), the Loc kman Hole bv | O wen & Morrisonl 
(12008 ), the COSMOS field by ISchinne rer et ah (120071), S tripe 
82 of the Sloan Digital Sky Sur vey bv I Hodge et al ns and 
the First-Look Surv ey region by iMorgantTetakl (l2004l i ~and by 
ICondon et alJ lUooT . The telescope used to conduct each of the 
surveys is indicated in the top right. The horizontal axis is the 5<r 
detection limit of the surveys in mjy and the vertical axis is the 
.S’foM factor for surveys in deg 2 mjy -2 (see section f2.II for more 
detail). The two dashed lines indicate the All Sky limit and the 
approximate envelope of existing surveys, heavily dependent on 
available telescope time. 


were then flagged using mirflag with int=1200 and further 
manually flagged with uvflag. The mean integration time 
per pointing after flagging is « 15 hours in CDFS and « 9 
hours in ELAIS-S1. 


2.3 Imaging 

2.3.1 Single Pointing Image Processing 

We combined all data for each pointing into one file using 
uveat and then each pointing was imaged separately. Since 
the observations cover a large fractional bandwidth, multi- 
freq uency synthesis was use d. 

ISault fe Conwavl (1 19991 ) describe the process of multi¬ 
frequency synthesis and how to successfully CLEAN an im¬ 
age. When creating an image over a wide bandwidth, two 
dirty beam images must be used, as the spectrum of a radio 


sources is given by 

I(v) = I(v 0 ) + al(v 0 )- —— , (2) 

vo 

where I is the flux density at frequency v, vp is the refer¬ 
ence frequency and a is the spectral index dSault fc Conwavl 
Il999l) . The two dirty beams are the synthesised dirty beam 
and the spectral dirty beam. The dirty image is then repre¬ 
sented by 

Id( l, m) =1(1, m ) * Bo(l, m) + ( a(l , m)I(l, m )) * Bi(l, m ) , 

(3) 

where Bo is the synthesised dirty beam, B 1 is the spec- 
tral dirty beam and ( l,m ) are the directional cosines 
dSault fe Conwavl Il999l ). The MIRIAD task invert provides 
the option of creating the two dirty beams. These two dirty 
beams are used in the multi-frequency CLEAN task, mf clean, 
to create the CLEAN component / and al maps. 

Each pointing was CLEANcd and self-calibrated sepa¬ 
rately prior to mosaicking. A model for self-calibration was 
produced by imaging the full band with uniform weight¬ 
ing and a cell size of 1 arcsec. The band was then split up 
into four 125 MHz subbands and self-calibration was ap¬ 
plied to each subband separately in order to allow for varia¬ 
tions in the antenna gains with frequency. Three iterations 
of phase self-calibration were applied, progressively increas¬ 
ing the number of CLEAN components used to model the sky 
emission, and hence the total amount of flux included in the 
model. 

The final self-calibrated uv data were then divided into 
two subbands of 250 MHz, imaged separately and adjusted 
to the reference pointing. Splitting up the data into two sub¬ 
bands serves to improve the accuracy of the primary beam 
correction. In order to obtain nearly identical synthesised 
beams for a ll pointings, individual ‘robust’ weighting factors 
dBriggsll 19951 ) were assigned. Pointings in the lower subband 
typically have a robust weighting factor in the range 0.0-0.5 
and in the higher subband 0.4-1.0. 

Using robust weighting factors in these ranges was 
found to minimise the rms noise in the final mosaic. Using 
higher robust weighting factors - a robust weighting factor 
of 2 or more corresponds very closely to natural weighting - 
was found to degrade the image sensitivity, despite the de¬ 
crease in the theoretical noise level. This probably arises, in 
large part, from the degradation of the sidelobe levels and 
beam-shape, which renders the task of CLEANing the image 
and removing artifacts harder. 

The larger beam size obtained when using natura l 
wei ghting also increases the confusion noise (ICondon et al.l 
120121) . The final beam size in ELAIS-S1 is 12.2 by 7.6 arc- 
sec and in CDFS is 16.3 by 6.8 arcsec (see Section 12.411 . 
At a frequency of 1.4 GHz and a resolution of 9.6 arcsec 
(the geometric mean of the major and minor beam axes in 
ELAIS-S1), the confusion noise cr c « 3.8 /iJy beam -1 . cr c is 
predicted to vary as 6> 10 / 3 , where 9 is the beam size. Using 
a robust weighting factor of 2 for all pointings in ELAIS-S1 
results in a final beam size of 17.9 by 11.8 arcsec, for which 
a c s=s 14.9 p Jy beam -1 . 

When imaging the self-calibrated uv data, the cell size 
was set to 1.5 arcsec. We CLEANcd each pointing down to 5cr. 
A 2D elliptical Gaussian fitted to the central region of the 
dirty beam was used as the restoring beam. The mf s option 
in restor was used to write a second plane in the output 
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Figure 2. Map of a 1.4 Jy source (NVSS J032836-284147) in 
CDFS, before (top) and after (bottom) amplitude self-calibration. 


image containing the al model convolved with the Gaussian 
beam. The al plane was subsequently used to perform wide¬ 
band primary beam correction in linmos (see Section m . 
The pointing reference for ELAIS-S1 is a = 00 h 33 m 50.8 s , 
5 = -43°44'57.4" and for CDFS is a = 03 h 30 m 16.3 s , 
5 = -28°05'12.4". 

2.3.2 Image Artifacts 

Radial spokes are present around bright sources in both 
fields. We found that these can be effectively removed us¬ 
ing amplitude self-calibration. Fig.[3]shows the improvement 
in dynamic range after amplitude self- calibration around a 
1.4 Jy source (NVSS J032836-284147; ICondon et al.lll998l l 
in CDFS; the dynamic range (as measured by dividing the 
peak flux density of the source by the rms in the vicinity 
of the source) is « 70 per cent higher after amplitude self¬ 
calibration. 

We tried to apply amplitude self-calibration to the AT¬ 
LAS pointings using components with flux density greater 
than 5a to model the sky brightness distribution. This 
caused the flux densities of sources below approximately 
1 mjy to be biased low. The magnitude of the bias was 


found to increase for sources with decreasing flux density 
below 1 mjy, reaching ~ 10 per cent close to the 5<r de¬ 
tection limit. We suspect that a considerable amount of 
flux was missing in the model due to the flat slope of the 
counts at sub-mjy level s, resulting in the obs erved bias in 
the source flux densities. IWieringa et al] (1 19921 ') investigated 
self-calibration methods in use at radio synthesis arrays and 
found that self-calibration can bias the gains if a signifi¬ 
cant amount of flux is missing in the model, particularly in 
cases where the number of antennas is low. Any sources not 
fully contained in the model tend to get absorbed into the 
gains and are then reduced in amplitude in the image after 
self-calibration. Amplitude self-calibration was therefore not 
used when producing the final images. 

Mild artifacts in the form of concentric rings are also 
visible around the strong source shown in Fig. [5] These lim¬ 
itations are related to difficulties in modelling and calibrat¬ 
ing time-dependent effects simultaneously using standard 
self-calibration techniques, particularly where the source is 
bright and partially (or fully) resolved. This is primarily due 
to the east-west nature of ATCA array configurations (i.e. 
at any point in ti me only a slice in the uv plane is sampled). 
iLenc et alj (l2009l f found that, in the pointing closest to the 
strong source, these artifacts could be removed by modelling 
the source in the uv plane with a combinat ion of Gaussian s 
and point-like components, using Difmap dShepherdlll997l f . 
We did not apply this technique here because it was difficult 
to automate and it could not handle wideband data. 


2.3.3 Addition of Previous ATLAS Data 

Previous ATLAS DR2 observations (ATC A Project IDs 
Cl035 and Cl241) from lHales et al] (l2014a|j were combined 
with our CABB observations to maximise the sensitivity. 
Information on these observations including observing dates 
and ATCA array configurations is given in Table I All The 
data were combined in the image plane for reasons related to 
bandwidth smearing, and described in detail in Section T3.21 
These data consist of two 128 MHz bands centred at 1344 
and 1432 MHz each containing 16 channels of 8 MHz in size. 
Both bands were imaged together using a similar procedure 
to the CABB data. The robust weighting factors for the 
pointings were chosen to yield a similar beam size to that of 
the CABB data. Pointings typically have a robust weighting 
factor in the range —0.1 —1-0.1. 

Addition of the pre-CABB data results in an improve¬ 
ment in sensitivity of ss 20 per cent and also improves the 
dynamic range thanks to the increased uv coverage. 


2.4 Mosaics 

Each pointing was convolved with a Gaussian to obtain an 
identical synthesised beam across each of the two fields sep¬ 
arately. CDFS has a beam size of 16.3 by 6.8 arcsec, with 
position angle 2°, and ELAIS-S1 a beam size of 12.2 by 
7.6 arcsec, with position angle —11°. A source-free estimate 
of the noise in each pointing image was obtained as follows: 
an initial estimate of the noise was obtained by taking the 
rms within the primary beam full-width at half-maximum 
(FWHM). To avoid the noise estimate from being affected 
by real source emission, all pixels outside the range ±3er 
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Figure 3. The ATLAS CDFS total intensit y mosaic with the linear greyscale set to the range —100 to +300/Dy beam -1 . The image 
projection used is North-celestial-pole (NCP; iGroise mu, a projection onto a plane perpendicular to the North Celestial Pole, which 
is a special case of the orthographic (SIN) projection, often used for east-west radio interferometers. The solid black contour indicates 
the component catalogue boundary (3.6 deg 2 ) of the mosaic defined by: (1) local rms noise + 100 //.Jy beam -1 ; (2) sensitivity loss due 
to bandwidth smearing < 20 per cent; and (3 ) mosaicked primary beam response p 40 per cent. The pattern of pointings on the sky is 
identical to that shown in lNorris et al.l ll2006i f. 


were then flagged and the rms was re-evaluated. This pro¬ 
cess was repeated a number of times until the noise was 
found to decrease by less than 10 per cent after removing 
pixel outliers. 

Finally, using linmos, the pointings for each field were 
corrected for the primary beam and mosaicked together, 
weighting them by their respective rms noise values; Gaus¬ 
sian primary beam fits for the new 16 cm CABB receivers 
across the entire frequency range (1.1 to 3.1 GHz j3 were used 
and the primary beam response was averaged over the sub¬ 
band as described below. Fig. [3] shows the resulting CDFS 
total intensity (Stokes I) image and Fig.[4]shows the result¬ 
ing ELAIS-S1 total intensity image. 

Due to the wideband nature of the observations there is 
a discrepancy between the integrated flux density over the 
band (as returned by MIRIAD) and the monochromatic flux 
density at the central frequency for any source with a non¬ 
zero spectral index. This is essentially due to most sources 


being best described by a power-law slope across the band 
rather than a simple linear slope. The integrated flux den¬ 
sity for a source with a power-law rather than a linear slope 
across the band is always going to exceed the monochro¬ 
matic flux density at the central frequency. For small ^ 
this effect is small, but increases for wider bands and for a 
source with an increasingly non-flat spectrum. In our case, 
for a = —0.75, typical of radio sources whose emission is 
dominated by optically thin synchrotron radiation, a source 
would have its monochromatic flux density over-estimated 
by less than 0.5 per cent in each subband, well within the 
absolute calibration errors. 

However, this issue gets more complicated off-axis as 
the attenuation of the primary beam drops off more rapidly 
at higher frequencies effectively reducing the sensitivity at 
higher frequencies. For narrow-band observations with a well 
known primary beam pattern (which varies with frequency), 
using the primary beam correction factor at the central fre¬ 
quency is sufficient. For wide-band observations the primary 


http://www.narrabri.atnf.csiro.au/people/ste616/beamshapes/bcamsha$)m.hSinEehtioh at this frequency is not truly representative. 
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Right Ascension (J2000) 

Figure 4. The ATLAS ELAIS-S1 total intensity mosaic with the linear greyscale set to the range —100 to +300/rjy beam -1 . The image 
projection used is NCP. The solid black contour indicates the component catalogue boundary (2.7 deg 2 ) of the mosaic defined by: (1) 
rms noise + 100 //.Jy beam -1 ; (2) sensitivity loss due to bandwidth smearing < 20 per cent; and ( 3) mo saicked primary beam response 
+ 40 per cent. The pattern of pointings on the sky is identical to that shown in lMiddclberg et alj J2008l j . 


To improve the accuracy of the primary beam cor¬ 
rection, we used the implementation of wideband primary 
beam correction in linmos (option bw): the image and pri¬ 
mary beam were evaluated at 10 frequencies, the image 
using the al plane and the beam using the beam fits 
across the frequency range. All these images and beams 
were then used in the standard mos aic equation (see e.g. 
ISault. Stavelev-Smith &: Brouwlll996l ) to produce the final 
mosaicked image. 

3 CORRECTION OF FINAL MOSAICS 
3.1 CLEAN Bias 

CLEAN bias dCondon et al.|[l998i ) is an effect in deconvolution 
which redistributes flux from point sources to noise peaks in 
the image, thereby resulting in a systematic underestimation 
of the flux densities of real sources. As the amount of flux 
which is taken away from a real source is independent of 
its flux density, the fractional error this causes is largest for 
faint sources. The magnitude of the effect will depend on 
the uv coverage and to what flux density level the images 
are CLEANcd. 

In order to analyse the CLEAN bias on the flux density of 


Table 2. Number of simulated sources added to the uv data of 
one ELAIS-S1 pointing at a given SNR to analyse the CLEAN 
bias. 


N 

SNR 

40 

5 

15 

6 

15 

7 

15 

8 

15 

9 

10 

10 

10 

12 

5 

16 

3 

20 

2 

30 

1 

50 

1 

100 


sources in our ATLAS mosaics, we followed a sim ilar proce¬ 
dure to that outlined in lMiddelberg et al.l (120081) . We added 
132 point sources at random positions to the uv data of one 
ELAIS-S1 pointing, except that a simulated source could not 
lie closer than 1 aremin from a real source (> 5<r) and any 
other simulated source. The numbers of sources added with 
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Figure 5. The mean normalised flux density of simulated sources 
as a function of the SNR. The lines illustrate the CLEAN bias when 
CLEANing to 5 a (black line and circles) and 2.5a (red line and 
squares). 


different signal-to-noise ratios (SNRs) are listed in Table [2] 
We then imaged the data in the same way as the final image 
(see Section 12.31) . extracted the flux density of each simu¬ 
lated source and divided by the injected flux density. This 
simulation was repeated 30 times to improve the statistical 
sampling. 

The mean normalised flux density of the simulated 
sources as a function of SNR is shown in Fig. [5] The flux den¬ 
sities of the simulated sources were measured by extracting 
the map values at the simulated source positions, interpo¬ 
lated between pixels. This was done by calculating the map 
values on a successively finer grid (up to 128 times finer), by 
repeated co nvolution with a Gaussian-graded sine function 
£r ee 111990). At low SNR, the thermal noise will introduce a 
slight shift in the position of a point source, resulting in an 
increased likelihood for the peak to lie on top of a positive 
noise fluctuation. This, in turn, causes the peak flux density 
of a point source to be biased slightly high; this is known as 
the peak flux density bias. We chose to measure the flux den¬ 
sities of the simulated sources at their true positions, rather 
than measure their peak flux densities, to disentangle the 
peak flux density bias from the CLEAN bias. The peak flux 
density bias is dealt with separately by the BLOBCAT source 
finder in Section \n\ 

The black circles in Fig. [5] show the CLEAN bias when 
CLEANing the images to 5a, which typically corresponds to 
5000 iterations (the final images were CLEANed to this level). 
The simulations indicate that a source with an SNR of 5 
would have its flux density underestimated by less than 1 per 
cent as a result of the CLEAN bias, well within the calibration 
error of 5 per cent (see Section 14.61) . We therefore chose not 
to correct the sources’ flux densities for this effect. 

We repeated the simulations CLEANing the images to 
2.5a (this typically corresponds to 65000 iterations) to see 
how the CLEAN bias is affected by the flux density level to 
which the images are CLEANed. In this case, the CLEAN bias 
is much more significant as shown by the blue squares in 



Speak/S (%) 

Figure 6. The fraction of the CDFS (solid line) and ELAIS-S1 
(dashed line) mosaics at or below a given bandwidth smearing 
level. The total area of the mosaic (O tota i) is defined where the 
rms noise level is less than 100 //Jy beam 1 . 


Fig. [5] the CLEAN bias is ~5 per cent for sources with an 
SNR of 5. 

3.2 Bandwidth Smearing 

Bandwidth smearing (chromatic aberration) is characterised 
by a smearing of the visibilities of a point source in the uv 
plane due to the range of frequency being sampled over a 
given bandwidth. This results in a decrease in the peak flux 
density of the source (ICondon et al.1 Il998l ) . The effect on 
sources in an image is to introduce smearing in the radial 
direction from the pointing centre, so although the peak 
flux density is reduced, the integrated flux density remains 
the same. Bandwidth smearing can cause the peak flux den¬ 
sity to fall below the sensitivity threshold, thus rendering it 
unrecoverable. The magnitude of the effect will depend on 
the source distance from the pointing centre and the frac¬ 
tional channel bandwidth. The effect of bandwidth smearing 
is small for CABB data (with 1 MHz channels) but signifi¬ 
cant for pre-CABB data (with 8 MHz channels) and is the 
primary c ause of the amp litude calibration discrepancy re¬ 
ported bv lNorris et al.l (120061 1 . 

In order to analyse the bandwidth smearing in our 
ATLAS mo saics, we f ollo wed a similar procedure to that 
outlined in lHales et al] (l2014al ) . We modelled the amount 
of bandwid th smearing to produce simulated images from 
lHales et al.l : 

•Speak _ l\ 21n2 ( A^eff d \1“ 5 / .\ 

s ~ [ 3 V F SproKcJJ ’ 1 j 

where S pea k is the peak flux density of the source, S the in¬ 
tegrated flux density of the source, Az/eff the effective chan¬ 
nel bandwidth, v the central frequency, d the radial dis¬ 
tance from the phase centre and _B pr0 j (£) the projected beam 
FWHM for an elliptical beam. The projected beam FWHM 



































ATLAS Third Data Release 9 



Figure 7. The fraction of the CDFS (solid line) and ELAIS-S1 
(dashed line) mosaics at or below a given rms noise level. The 
total area of the mosaic (f^total) is defined where the rms noise 
level is less than 100 jt/Jy beam -1 . 


is given by [Hales et al.l : 
-Hproj(C) — 


-Hmaj -Hm in 


\/[-BmajSill(£ - Ip )] 2 + [_B min COs(C — Ip )] 2 


(5) 


where £ is the position angle (East of North) of a source with 
respect to the phase centre, B ma j and B m i n are the major 
and minor axes FWHM of the elliptical beam respectively, 
and ip is the beam position angle measured East of North. 

All of these simulated images were mosaicked together 
following the method outlined in Section [2.41 The effect of 
bandwidth smearing on the CDFS and ELAIS-S1 mosaics 
is shown in Fig. [6] The bandwidth smearing for both mo¬ 
saics is less than 10 per cent over the great majority of the 
observational area. Sources’ peak flux densities were cor¬ 
rected for bandwidth smearing, using our maps quantifying 
the sensitivity loss due to bandwidth smearing, as described 

in Section 521 


4 1.4 GHz COMPONENT CATALOGUE 

4.1 Background Noise Map 

In order to extract sources from the mosaics we require an 
estimate of the background noise level. The background rms 
noise of each mosaic was c alculated for each p ixel follow¬ 
ing the method described bv lFranzen et all (l2014ll . For each 
pixel, the noise was taken as the rms inside a box of size 
21 + 1 centred on the pixel; we chose the half-width of the 
box, l, to be 20 times the beam for both mosaics to m inimise 
the number of spurious detections (I Huynh et al.ll2012h . How¬ 
ever, in order to avoid the noise estimate from being affected 
by real source emission, points were clipped iteratively until 
convergence at ±3er was reached. Fig. [7] shows the fraction 
of the ATLAS sky at or below a given rms noise level, where 
the ATLAS sky is defined where the rms noise is less than 


Table 3. Statistics for the rms noise distribution across the AT¬ 
LAS sky (region where the rms noise is less than 100 //.Jy beam -1 ) 


for CDFS and ELAIS-S1. 



CDFS 

ELAIS-S1 

Mean (/xjy beam -1 ) 

30 

35 

Median (fiJy beam -1 ) 

23 

24 

Mode (/xJy beam -1 ) 

14 

17 


100 /rJy beam 1 . Tabic [4] gives statistics for the rms noise 
distribution across the ATLAS sky for the two fields. 


4.2 Source Extraction 


We limited the area of the two ATLAS fields for source ex¬ 
traction defined by the union of the following criteria: 

(1) rms noise ^ 100 pJy beam -1 . 

(2) sensitivity loss due to bandwidth smearing < 20 per 
cent. 

(3) mosaicked primary beam response > 40 per cent of 
the peak response. 


Fig.[3]shows the defined area covering 3.6 deg 2 of the CDFS 
mosaic and Fig.[I]shows the defined area covering 2.7 deg 2 of 
the ELAIS-S1 mosaic. The resulting area is primarily defined 
by the primary beam response. 

We extracted source co mponents from the two AT¬ 
LAS fields using B LOBCAT jHales et al.1 l2012lh AEGEAN 
dHancock et al.l l2012h, pyBDsT and the AMI S 0URCE_FIND 
software ( AMI Consortium: Franzen et al.ll201l|j . When us¬ 
ing BLOBCAT, AEGEAN and pyBDSM, we searched both mosaics 
down to an SNR ^ 4 in order to include all sources with 
fitted peak SNRs ^ 5. Sources with fitted peak SNRs < 5 
were removed manually from the output component lists. 
When running SOURCE_FIND on the two mosaics, a detec¬ 
tion threshold of 5cr was used since the detection threshold 
is automatically lowered to ensure that all sources with fit¬ 
ted peak SNRs ^ 5 are included. We compared the source 
finders by locating the corresponding counterparts to within 
20 arcsec. The difference between the four source finders 
for our work is negligible. However, we chose BLOBCAT over 
the other source finders as BLOBCAT takes into consideration 
bandwidth smearing and peak bias. A more detailed anal¬ 
ysis of the source finders is beyond the scope of this paper 
but is discussed by Hopkins et al. (2014), submitted. 

Components that corresponded to artifacts around 
bright sources were removed from the BLOBCAT catalogues 
for the two fields, resulting in a total of 2861 components 
for CDFS and 1964 for ELAIS-S1. 


BLOBCAT assumes that isolated components have Gaus¬ 
sian morphology in order to catalogue their properties. This 
assumption may not always be suitable for complex com¬ 
ponents. We identified complex blobs as having A p i x ^ 300 
and Rest ^ 1-4, where N p i x is the number of flooded pixels 
comprising the blob and R es t is the size estimate of the blob, 
in units of the sky area covered by an unresolved Gaussian 
blob with the same peak flux density. Gaussian fitting was 


3 http://dl.dropboxusercontent.com/u/1948170/html/index.html 
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pe rformed for each of these complex blobs, as recommended 
bv lHales et al.l J2012I) , with the CASA task imfit. Initially, 
two Gaussians were fitted simultaneously and the residuals 
were examined. If the peak of the residuals was less than 5a, 
then the fit was considered to be acceptable and the original 
BLOBCAT catalogue entry was replaced by the imfit entry for 
each individual Gaussian component. Otherwise, the num¬ 
ber of Gaussians included in the fitting was increased by 
one and the residuals were re-evaluated. This process was 
repeated until the peak residual was found to be less than 
5 a or 8 Gaussians had been fitted. 

Following this procedure, Gaussian fitting was per¬ 
formed for a total of 97 and 72 BLOBCAT components identi¬ 
fied as complex in CDFS and ELAIS-S1, respectively. We 
then merged together those Gaussian components which 
were separated by less than half the beam size. Our final 
component catalogue for CDFS contains 3034 components 
and that for ELAIS-S1 contains 2084 components. 


4.3 Deconvolution 

A measure of the degree to which a radio source is resolved 
is given by the ratio of the integrated flux density, S, to the 
peak flux density, Speak: 


Speak 


^min^max 

B m aj Bmin 


(6) 


where # m i n and # max are the source fitted FWHM axes. 

In Fig. 0 we plot S/S pea k as a function of SNR for 
CDFS in the left panel and for ELAIS-S1 in the right panel. 
The peak flux density has been corrected for bandwidth 
smearing and peak bias. The distribution of components is 
skewed to higher flux density ratios at higher SNR, a re¬ 
sult of extended components. The occurrence of components 
with S/Speak < 1 are from noise fluctuations and calibration 
errors. For point sources, assuming that as peak and as are 
independent, R = ln(S'/S pe ak) has a Gaussian distribution 
centred on zero whose rms is given by 



To detect an extended source at the 2 a level, we require 
R > 2<tr or 



The probability of falsely classifying a point source as ex¬ 
tended at the 2a level is 2.3 per cent. 

We used Equation 0 to separate point-like from ex¬ 
tended sources. as peak is taken as the sum in quadrature of 
the calibration er ror, the pixellation uncertainty (see Appen¬ 
dices A and B of lHales et al]|2012h and the local rms noise; 
the calibration error is set to 0.05 Speak (see Section [4.6.11) 
and the pixellation uncertainty to 0.01 Speak- as is taken as 
the sum in quadrature of the calibration error and the local 
rms noise; the calibration error is set to 0.05 S. At high SNR, 
where as/S 0.05 and as peak /Speak « 0.05, S/S pe ak must 
be greater than 1.15 for a source to be classified as extended. 
The number of components classified as extended is 514 (17 
per cent) in CDFS and 392 (19 per cent) in ELAIS-S1. 


The deconvolved angul ar size for extended sources is 
given bv lHales et al.l d2014al l: 


0 — \/ ^min^max B m aj Rn 


(9) 


BLOBCAT does not directly measure (9 m i n or 0 ma x. By substi¬ 
tuting for # m in#max using Equation 0 



( 10 ) 


Following standard error propagation, the error on 0 is given 
by: 


a e = 


B m aj Rmin 


•Speak \ 4(S/Speak ~ 1) 


f ^Speak 

+ 

\ <^peak 

\ s) 


( 11 ) 


4.4 Spectral Indices 


The wide bandwidth of the ATLAS observations provides 
enough information to calculate the spectral indices over 
the mosaics. We have measured the spectral indices of com¬ 
ponents where the SNR was high enough to obtain a reliable 
measurement exclusively using the CABB data. 

We created two separate mosaics of each held, one us¬ 
ing the lower CABB subband data centred at v\ — 1.40 GHz 
and the other using the higher CABB subband data centred 
at V 2 = 1.71 GHz, ensuring that both mosaics had the same 
resolution. We note that flagged frequency channels were 
accounted for in MIRIAD when calculating the central fre¬ 
quencies of the two subbands. 

To measure the spectral indices of isolated components, 
we ran BLOBCAT on the two CABB subband mosaics using a 
5a detection limit. If a component was classified as extended 
in the final catalogue, we used its integrated flux density at 
ui and U 2 as measured by BLOBCAT to calculate its spectral 
index, otherwise we used its peak flux density. We did not 
measure the spectral indices of components detected below 
5a at u\ or ■ 

We derived the spectral indices of complex compo¬ 
nents by measuring their integrated flux densities in the 
two CABB subband images. The integrated flux density of 
a complex component in each subband image was measured 
by summing the pixel values within its integration area (as 
measured in the raw BLOBCAT catalogue described in Sec- 
tion l4.2l) . and dividing by the number of pixels per beam; the 
advantage of this method for measuring the spectral index 
is that it ensures that the flux density within the exact same 
region of the image is compared at each frequency, which is 
important for components spanning several beam areas. We 
did not attempt to measure the spectral indices of Gaussian 
components which were fitted to complex components. 

The error on the spectral index was taken as 


= \J <^,th + <^,cal , (12) 

where a^th is the uncertainty on a due to the thermal noise 
and a a , C ai is the calibration error on a. The noiselike uncer¬ 
tainty on a is given by 

T 


r 




/ 

ct ("i) 
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[sM\ 


|so 2 )J 
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(13) 
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Figure 8. S/S pegi k as a function of SNR for all components in CDFS (left) and ELAIS-S1 (right). The peak flux density values have 
been corrected for peak bias and bandwidth smearing as measured in BL0BCAT. Components which are classified as point-like/extended 
are shown in black/red. 



10 100 1000 10000 
SNR 



where S(ui) and S(v 2 ) are the flux densities at V\ and V 2 , 
a(u 1 ) and 17 ( 1 / 2 ) are the rms local noise at v\ and z/ 2 , and 
P = \xi(y 2 /v\). For complex components spanning multiple 
beam areas, we multiplied a(v 1 ) and ( 7 ( 1 / 2 ) by the square 
root of the integration area in units of the synthesised beam. 

We have identified two main types of calibration errors 
contributing to er aiCa i: uncorrelated errors in the primary 
beam model at and 1 / 2 , and uncorrelated errors in the 
antenna gain calibration at u\ and V 2 ■ Primary beam model 
errors are expected to have the largest contribution close to 
the edges of the mosaics. As mentioned in Section l4~2l source 
extraction was limited to the region where the mosaicked 
primary beam response is above 40 per cent. The 40 per cent 
power level at V 2 corresponds to a distance from the pointing 
centre, d, of approximately 16 arcmin. In Section [4.4.II we 
show that primary beam model errors cause a to flatten by 
no more than ss 0.1 at d = 20 arcmin. We conclude that, 
for a component located at any position within the mosaics, 
primary beam model errors introduce a spectral index error 
of at most rb 0.1. 

Given the very high degree of correlation (<; 95 per 
cent) of gain calibration errors in the two CABB subbands, 
we have established that their contribution to cr Q , C ai is much 
less than 0.1. Since primary beam model errors are likely to 
be the dominant contribution to ff Q , C ai, we set cr Q|Ca i = 0.1. 

In total, we measured spectral indices for 1756 isolated 
components and 163 complex components. Of the 1919 spec¬ 
tral index measurements, 344 have a a less than 0.2. We have 
included a column containing the spectral index in the final 
component catalogue for each field (see Table [5]|. The spec¬ 
tral index properties of the ATLAS DR3 components are 
discussed further in Paper II (Banfield et al., in prepara¬ 
tion). 

4-4-1 Verifying off-axis spectral index measurements 

In 2010, measurements of the primary beam for the new 
ATCA 16 cm receivers were carried out. The radial profile of 
the beam was measured along eight cuts at 45° intervals, at 
14 regularly-spaced frequencies across the band. The shape 
of the beam out to the first null was found to be close to 


circularly symmetric across the entire frequency range. A 
Gaussian fit to the beam was made at each frequency out to 
a level of ~ 20% of the beam, using data from all angles. The 
product of the primary beam FWHM and the frequency was 
found to increase by 7 per cent between 1.3 and 1.8 GHz. 
This is probably due to a slight defocussing between these 
two frequencies. 

As mentioned in Section El we used these Gaussian 
primary beam fits for primary beam correction. A small er¬ 
ror in the primary beam model may introduce a bias in the 
measured spectral indices away from the pointing centre. We 
therefore tested whether there was any systematic change in 
the spectral index with distance from the pointing centre. 
We applied the primary beam correction to each pointing in 
CDFS and ELAIS-S1 and restored them to a common res¬ 
olution. We identified sources in the raw BL0BCAT catalogue 
for the relevant field which were located within 20 arcmin 
from the pointing centre and with SNR > 200. We measured 
the spectral indices of these sources by comparing the pixel 
values at their peak positions in the two CABB subband 
images. The top panel of Fig. [9] shows the measured spec¬ 
tral index as a function of the distance from the pointing 
centre, combining the results from all pointings. Given the 
SNR. cut applied, the errors on the spectral indices are less 
than approximately 0.1. We measured the median spectral 
index in five distance bins (0-4, 4-8, 8-12, 12-16 and 16- 
20 arcmin) and fit a second order polymonial to the median 
data points. The median data points are shown in Table 3] 
There is no statistically significant change in the spectral in¬ 
dex with distance from the pointing centre out to 20 arcmin, 
which corresponds to a power point of rb 0.25 at 1.71 GHz. 
At d = 20 arcmin, we can be confident that errors in the pri¬ 
mary beam model cause a flattening in the spectral index of 
no more than rb 0.1. 

To illustrate how sensitive the spectral index is to er¬ 
rors in the primary beam model, the bottom panel of Fig. [9] 
shows how the spectral index varies with distance from 
the pointing ce n tre us ing an older primary beam model by 
IWierinaa et alj (1 1992h . where the beam FWHM varies as 
i/ _1 between 1.3 and 1.8 GHz. In this case, there is a clear 
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Figure 9. Spectral index between 1.40 and 1.71 GHz as a func¬ 
tion of distance from the pointing centre for bright sources present 
in the CDFS and ELAIS-S1 pointings. The red points show me¬ 
dian values in five distance bins, which are tabulated in Table 21 
The blue curve shows a quadratic fit to the median data points. 
Top: results using Gaussian primary beam fits for the new 16 cm 
CABB receivers between 1.1 and 3.1 GHz. Bott om: results ob¬ 
tained using an older primary beam model by IWieringa et al.l 
(119921) . where the beam FWHM varies as v 1 between 1.3 and 
1.8 GHz. 


bias in the spectral index away from the pointing centre; the 
bias is « 0.25 at the half power point at 1.71 GHz. 


4.5 Frequency Coverage 

The combination of multiple epochs of ATLAS observa¬ 
tions, the wider frequency coverage of the DR3 data and 
the amount of flagging resulted in the effective observing 
frequency changing slightly across the mosaics. We used the 
task linmos to produce a mosaic giving the effective fre¬ 
quency across each field. In the final component catalogue, 
we have included the frequency at which the source’s flux 
density is measured. The frequency typically lies in the range 
1.45-1.50 GHz. 

We have also provided a column with the source’s 
flux density at 1.4 GHz. This was derived using ai;Io if 
cr Q i.7i Sj 0.2, which is typically the case for sources with 
Si .4 > 2 mjy. In their stu dy of the sub-mjy radio popula¬ 
tion in the Lockman Hole, llbar et all (120091 ) found no sig¬ 
nificant change in the median value of as a function of 
Si. 4 ; aj;!? was found to be approximately -0.6 to -0.7 down 
to a flux density level of Si .4 > 100/xJy. For ATLAS sources 


Table 4. Relationship between spectral index and distance from 
the pointing centre for bright sources in individual CDFS and 
ELAIS-S1 pointings. The first column gives the distance range 
and the second column the number of sources. The third column 
gives the median spectral index obtained using Gaussian primary 
beam fits for the new 16 cm CABB receivers between 1.1 and 
3.1 GHz. The numbers in brackets are median spectral indices 
obtain ed using an older primary beam model by I Wieringa et akl 
1 19921) . where the beam FWHM varies as ^ 1 between 1.3 and 
1.8 GHz. 


r 

N 

^median 

0-4 

17 

-0.82 ±0.11 (-0.79 ±0.10) 

4-8 

59 

-0.78 ± 0.05 (-0.72 ± 0.06) 

8-12 

105 

-0.72 ± 0.04 (-0.65 ± 0.04) 

12-16 

159 

-0.74 ± 0.04 (-0.53 ± 0.05) 

16-20 

184 

-0.70 ± 0.04 (-0.38 ± 0.04) 



Figure 10. Ratio of the ATLAS to VLA flux density as a function 
of the VLA flux density for 112 unresolved s ources in the eCDFS . 
The VLA flux densities were obtained from [Miller et ah! J2013I) . 
The dashed horizontal line indicates equal flux density values 
and the dashed vertical line indicates the flux density (0.15 mJy) 
above which the ATLAS flux densities are not considered to be 
affected by the Eddington bias. 


with cr^i.71 > 0.2 or with no measured spectral indices, we 

a 1.40 

therefore assumed that aq^o = —0.7 to obtain Si. 4 . 

4.6 Calibration accuracy checks 

I Miller et al.l (|2013l ) produced an image of the Extended 
Chandra Deep Field South (eCDFS) with the VLA at 
1.4 GHz. Their image covers an area of about a third of 
a square degree to a typical rms sensitivity of 7.4 /ijy bm -1 
and has a resolution of 2.8 by 1.6 arcsec. In order to as¬ 
sess the accuracy of the ATCA flux densities and posi¬ 
tions presented in this paper, we matched our catalogue 
wit h the mo re sensitive and higher resolution catalogue by 
iMiller et all 

4-6.1 Flux density calibration accuracy 

The flux den sity s cale for ATCA is consistent with that of 
I Baars et al.l (Il973l ). in use at th e VLA, at the 1 -2 per cent 
level over the range 1-10 GHz (iRevnoldsI 1 19941) . The VLA 
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Figure 11. RA and Dec offsets between the ATLAS and VLA 
positions of 112 unresolved sources in the eCDFS. Sources with 
ATCA SNRs ^ 20 are shown in red and the rest of the sources in 
black. 



is more likely to resolve out part of the source flux den¬ 
sity because of its smaller beam size. In order to minimise 
discrepancies between the ATCA and VLA flux densities re¬ 
sulting from the smaller VLA beam size, we only considered 
sources which were classified as point-like in both the ATCA 
and VLA catalogues. 

Fig. 1101 shows how the ATCA and VLA integrated flux 
densities compare for all 112 unresolved sources in common 
between the two surveys. There is generally excellent agree¬ 
ment between the flux densities. At the faint end, the ATCA 
flux densities tend to be systematically higher than the VLA 
flux dens ities. This is prob ably due, in part, to the Edding¬ 
ton bias (lEddingtonll 191 j l causing the ATCA flux densities 
to be biased high close to the survey detection limit; the rms 
noise in the eCDFS region of ATLAS is ~ 3 times higher 
than that in the VLA image. The flux density boosting re¬ 
sulting from Eddington bias only depends on t he SNR and 
the sou rce count slope dHogg fe Turned [l998lb Ildales et al~ 
(20Ma|) used the 1.4 GHz source count fit bv lHopkins et al. 
( 20031) to evaluate the flux density boosting for a 5<r source 
as a function of Si. 4 ; the bias was estimated to be « 10 
per cent at Si .4 = 0.1 mjy. This is not sufficient to explain 
the observed discrepancy between the ATCA and VLA flux 
densities at the faint end. Another potential cause of the 
discrepancy is missing extended flux density in the VLA 
image. 

For sources with Svla > 0.15 mjy, the ATCA flux den¬ 
sities do not appear to be affected by the Eddington bias. 
For these 74 sources, the ATCA flux densities are, on aver¬ 
age, consistent with the VLA flux densities at the « 2 per 
cent level: the mean value of Satca/SVla is 0.986 ± 0.013 
and the median value is 0.974 ± 0.016. 

At high SNRs, calibration errors will become the main 
source of uncertainty. Of the 112 ATCA-VLA sources, 29 


have ATCA SNRs ^ 20. For these 29 sources, the stan¬ 
dard deviation of Satca/SVla is 0.047. We therefore set 
the ATCA flux density calibration error to 5 per cent. This 
is a conservative estimate of the flux density calibration er¬ 
ror since it does not account for the scatter in Satca/Svla 
introduced by errors on the VLA flux densities. 


4-6.2 Positional accuracy 

We have estimated the calibration errors a a}CS i and as :C ai 
in RA and Dec by comparing the ATCA positions with the 
VLA positions. Fig. llll shows the RA offset, d a , and the Dec 
offset, d-s, for the 112 unresolved sources in common between 
the two surveys. Sources with ATCA SNRs ^ 20, for which 
calibration errors dominate the position uncertainties, are 
shown in red. 

The positional uncertainties in RA and Dec resulting 
from phase errors can be expressed as b Q a and b Qg, where 
b is a constant, and Q a /g is the projected resolution in the 
RA/Dec direction. For the 29 sources with ATCA SNRs 
^ 20, the standard deviation of d a is 0.13 arcsec and the 
standard deviation of ds is 0.36 arcsec, indicating that b « 
0 . 021 . 

To test for systematic errors, we calculated the mean 
offsets in RA and Dec. The mean value of d a is (—0.12 dt 
0.02) arcsec and the mean value of dg is (0.07 ± 0.07) arc¬ 
sec. To account for the systematic offset in RA, we in¬ 
creased <j a,cal by 0.12 arcsec. In CDFS, we set <r a ,cai = 
b 0q. + 0.12 arcsec = 0.26 arcsec. Since the mean value of 
dg is consistent with zero at the ler level, we assumed no 
systematic offset in Dec and set crj jC ai = b Qg = 0.33 arcsec. 

It was not possible to carry out a similar analysis in 
ELAIS-S1 because no high resolution radio data were avail¬ 
able to perform a comparison. Assuming no systematic er¬ 
rors in RA and Dec, and using the same value of b as for 
CDFS, cr ajCa i = b 0 Q = 0.16 arcsec and crg iCa \ = b Qg — 
0.25 arcsec. 

We set the total position errors cr Q and ag in RA and 
Dec to 

= \J oi, cal + °i, th (14a) 


o-g 



+ CT 


2 

<5,th ’ 


(14b) 


where cr aj th and og : th are the position errors in RA and Dec 
due to the significance of the component detection, blobcat 
(Hales et al. 2012) estimates cr a ,th and erg ,th as 


& Q,th 


1.4 A 


O a 


(15a) 


& <5,th 


where A is the SNR. 



(15b) 


4.7 Component Catalogue 

We have combined the source component catalogues for both 
the CDFS and ELAIS-S1 into one component source cata¬ 
logue for ATLAS. Table [5] fists the first 14 components in 
the catalogue while the full catalogue can be obtained from 
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the supplementary material. The columns of the component 
catalogue are defined as follows: 

Column (1) - component number. A prefix indicates 
the field in which the component lies: ‘CF and ‘EF stand for 
Stokes I component in CDFS and ELAIS-S1, respectively. 
Where multiple Gaussians were fitted to complex compo¬ 
nents, the component number is appended by the letter ‘C’ 
followed by the Gaussian number. 

Column (2) - component IAU name given by ATLAS3 
Jhhmmss.s-ddmmssC, where the letter ‘C’ stands for com¬ 
ponent in Stokes I. Norris et al, in preparation compiles 
the ATLAS DR3 source catalogue by combining the source 
components. 

Columns (3) and (4) intensity weighted centroid po¬ 
sition: RA (J2000), a, in hours:minutes:seconds, and Dec. 
(J2000), 5, in degrees:minutes:seconds. 

Columns (5) and (6) intensity weighted centroid po¬ 
sition: RA (J2000), a, and Dec. (J2000), S, in deg. 

Columns (7) and (8) - error on centroid position in 
RA, a a , and in Dec., as, in arcsec. These were derived as 
described in Section T4. 6.2 1 

Column (9) local rms noise level, criocai, in mjy 
beam -1 . 

Column (10) - bandwidth smearing correction, b (see 
Section IT21) . 

Column (11) ~ effective frequency, r'obs, in GHz, at 
which the source was measured in the mosaic (see Sec¬ 
tion 14.511 . 

Columns (12) and (13) - peak flux density, Speak, and 
associated error, <rs peak , at v 0 b B , in nrjy beam -1 . S pea k has 
been corrected for peak bias and bandwidth smearing in 
BLOBCAT. a s peak is derived as described in Section ^obs 
is provided in column (20). 

Columns (14) and (15) - integrated flux density, S, and 
associated error, as, at i^obs, in mjy. as is derived as de¬ 
scribed in Section [431 S is set to S pea k and as to as peak for 
sources classified as point-like. 

Columns (16) and (17) - peak flux density, S pe ak,i. 4 , 
and associated error, o\g peak 1 4 , at 1.4 GHz, in mjy beam -1 . 
•Speak, 1.4 was derived from Speak as described in Section 1431 

Columns (18) and (19) - integrated flux density, Si. 4 , 
and associated error, o\g 14 , at 1.4 GHz, in mjy. Si .4 was 
derived from S as described in Section S3] 

Columns (20) and (21) - deconvolved angular size, 0, 
and associated error, ae, in arcsec, for sources classified as 
extended. 0 and ae are set to —999 for point sources. 

Column (22) - component type, t\ point-like (P) or ex¬ 
tended (E) (see Section S3]). 

Columns (23) and (24) - spectral index between 1.40 
and 1.71 GHz, al'do, and associated error, cr a i.7i, as mea¬ 
sured in SectionSIU The spectral index is set to —999 if the 
SNR was too low for it to be measured. 

Column (25) - field identifier. 




Table 5. Catalogue entries for the first 14 components in the ATLAS 1.4 GHz component catalogue. The columns are defined in Section l4.7j The full catalogue can be found in the 
supplementary material using the same column format. 


Component number 

(i) 

IAU Name 

(2) 

a 

(3) 

6 

(4) 

a 

(deg) 

(5) 

S 

(deg) 

(6) 

(T a 

(arcsec) 

(7) 

(arcsec) 

(8) 

^local 
(mjy bm —1 ) 

(9) 

b 

(10) 

^obs 

(GHz) 

(ii) 

*Speak 
(mjy bm -1 ) 
(12) 

EI0001 

ATLAS3 J002925.7-440256C 

00:29:25.73 

-44:02:56.6 

7.357191 

-44.049053 

0.44 

0.69 

0.052 

1.08 

1.439 

0.751 

EI0002 

ATLAS3 J002926.7-440016C 

00:29:26.79 

-44:00:16.5 

7.361610 

-44.004578 

1.07 

1.65 

0.051 

1.08 

1.440 

0.288 

EI0003 

ATLAS3 J002933.7-440118C 

00:29:33.70 

-44:01:18.4 

7.390433 

-44.021768 

0.79 

1.22 

0.046 

1.08 

1.447 

0.349 

EI0004 

ATLAS3 J002938.1-432946C 

00:29:38.17 

-43:29:46.2 

7.409025 

-43.496162 

0.34 

0.52 

0.045 

1.09 

1.446 

0.904 

EI0005 

ATLAS3 J002938.9-440031C 

00:29:38.91 

-44:00:31.2 

7.412115 

-44.008663 

0.70 

1.08 

0.043 

1.07 

1.452 

0.367 

EI0006 

ATLAS3 J002940.1-440308C 

00:29:40.19 

-44:03:08.7 

7.417454 

-44.052408 

0.60 

0.92 

0.042 

1.07 

1.452 

0.434 

EI0007 

ATLAS3 J002941.3-435334C 

00:29:41.34 

-43:53:34.0 

7.422268 

-43.892779 

0.61 

0.95 

0.043 

1.08 

1.440 

0.430 

EI0008 

ATLAS3 J002941.8-440714C 

00:29:41.80 

-44:07:14.1 

7.424153 

-44.120575 

0.99 

1.53 

0.042 

1.08 

1.447 

0.252 

EI0009 

ATLAS3 J002943.1-440813C 

00:29:43.11 

-44:08:13.3 

7.429629 

-44.137028 

0.77 

1.19 

0.042 

1.08 

1.446 

0.333 

EI0010 

ATLAS3 J002944.3-433629C 

00:29:44.30 

-43:36:29.0 

7.434567 

-43.608068 

0.63 

0.98 

0.049 

1.09 

1.435 

0.480 

EI0011C1 

ATLAS3 J002945.3-432148C 

00:29:45.34 

-43:21:48.0 

7.438908 

-43.363342 

0.16 

0.26 

0.042 

1.08 

1.444 

6.506 

EI0011C2 

ATLAS3 J002946.1-432148C 

00:29:46.17 

-43:21:48.7 

7.442358 

-43.363537 

0.16 

0.26 

0.042 

1.08 

1.445 

7.060 

EI0012 

ATLAS3 J002946.3-440724C 

00:29:46.33 

-44:07:24.0 

7.443055 

-44.123333 

0.82 

1.27 

0.040 

1.07 

1.450 

0.289 

EI0013 

ATLAS3 J002948.6-435618C 

00:29:48.66 

-43:56:18.6 

7.452745 

-43.938512 

0.71 

1.10 

0.037 

1.07 

1.454 

0.315 


Table 5 - continued 


'’’Speak 

S 


‘Speak,].. 4 

^^peak,!.^ 

Si.4 

<T Sl.4 

© 

o-q 

t 

rv 1 - 71 

a 1.40 


Field 

(mjy bm -1 ) 

(mJy) 

(mjy) 

(mjy bm -1 ) 

(mjy bm -1 ) 

(mJy) 

( m Jy) 

(arcsec) 

(arcsec) 




(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

0.068 

0.751 

0.068 

0.766 

0.070 

0.766 

0.070 

-999.0 

-999.0 

P 

-0.83 

1.01 

ELAIS-S1 

0.057 

0.288 

0.057 

0.294 

0.059 

0.294 

0.059 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.052 

0.349 

0.052 

0.357 

0.054 

0.357 

0.054 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.067 

0.904 

0.067 

0.924 

0.069 

0.924 

0.069 

-999.0 

-999.0 

P 

-2.72 

0.93 

ELAIS-S1 

0.049 

0.367 

0.049 

0.377 

0.051 

0.377 

0.051 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.051 

0.434 

0.051 

0.445 

0.052 

0.445 

0.052 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.051 

0.430 

0.051 

0.438 

0.052 

0.438 

0.052 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.047 

0.252 

0.047 

0.258 

0.048 

0.258 

0.048 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.049 

0.333 

0.049 

0.340 

0.050 

0.340 

0.050 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.059 

0.480 

0.059 

0.489 

0.060 

0.489 

0.060 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.330 

7.983 

0.406 

6.703 

0.341 

8.224 

0.420 

4.6 

0.9 

E 

-0.96 

0.15 

ELAIS-S1 

0.357 

8.632 

0.438 

7.278 

0.370 

8.899 

0.454 

4.5 

0.9 

E 

-0.96 

0.15 

ELAIS-S1 

0.045 

0.289 

0.045 

0.296 

0.046 

0.296 

0.046 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 

0.043 

0.315 

0.043 

0.323 

0.044 

0.323 

0.044 

-999.0 

-999.0 

P 

-999.00 

-999.00 

ELAIS-S1 


Oi 


ATLAS Third Data Release 








16 T. M. 0. Franzen et al. 


5 CONCLUSION 

We present images and a component catalogue from a deep, 
wideband, radio continuum survey with ATCA. The third 
ATLAS data release combines observations taken between 
2002 and 2010 of an area coincident with the CDFS and 
ELAIS-S1, two of the best-studied regions of the sky at all 
wavelengths. The survey covers a total area of 6.3 deg 2 to a 
typical rms noise level of 15/rJy beam -1 at 1.4 GHz. Vari¬ 
ous array configurations were employed to maximise the uv 
coverage, resulting in a resolution of 16 by 7 arcsec in CDFS 
and 12 by 8 arcsec in ELAIS-S1. ATLAS is among the deep¬ 
est and widest radio surveys to date and is being used as 
a pilot survey for EMU, which will cover the whole South¬ 
ern Sky to approximately the same depth as ATLAS, at a 
similar resolution and frequency. 

In Paper II, we present the first results from the survey, 
including the deep 1.4-GHz source counts. Here, we have 
concentrated on developing techniques for producing and 
analysing the radio maps to enable maximum scientific re¬ 
turn from the survey. In particular, we have: 

(1) developed strategies to automatically flag data taken 
with the new CABB correlator between 1.3 and 1.8 GHz, 
and calibrate the data while accounting for frequency- 
dependent variations in the gains; 

(2) used a variety of techniques, such as multi-frequency 
CLEAN, self-calibration and wideband primary beam cor¬ 
rection, to produce images with high dynamic range and 
fidelity; 

(3) run the BLOBCAT source finder on our maps while ap¬ 
plying corrections for bandwidth smearing and peak bias to 
the source flux densities, resulting in a total of 5118 compo¬ 
nents above 5a in the two fields. 

(4) used the wide bandwidth of our observations to mea¬ 
sure the spectral indices across the CABB band (1.3- 
1.8 GHz) of ~ 2000 of the brightest components in the fields. 
Analysis of the spectral index results will be presented in 
Paper II. 
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Table Al. Observational information on the two ATLAS fields including observing dates, ATCA array configurations, and net integration 
times after calibration and flagging. 


Field Project 

Date 

Configuration 

Integration 

ID 



Time (h) 

CDFS C1035“’ c ’ d 

2002 Apr 4-7, 10, 12-13 

6A 

72.9 


2002 Aug 23-24, 27-29 

6C 

29.6 

Cl241 a ’ c ' d 

2004 Jan 7-8, 12 

6A 

23.9 


2004 Feb 3-5 

6B 

24.7 


2004 Jun 6, 8-12 

750D 

37.4 


2004 Nov 24-30 

6D 

50.4 


2004 Dec 28-30 

1.5D 

22.6 

C1241 c ’ d 

2005 Jan 7-8, 18-19, 23 

750B 

31.9 


2005 Apr 9-10 

6A 

18.5 


2005 Apr 14 

1.5A 

8.9 


2005 Apr 22; 2005 May 2 

750A 

15.0 


2005 Jun 1, 10 

EW367 

11.7 


2005 Jun 25-26 

6B 

18.1 


2005 Dec 6 

6A 

8.7 


2006 Mar 23-24, 27 

6C 

23.0 

C1967 d 

2009 Jun 19-21, Dec 19-31, 2010 Jan 1-3, Apr 1-6, 12, 14-18 

6A 

148.1 


2009 Aug 13, 19, 21-27 

6D 

76.6 


2009 Nov 10-18, 20-29 

6B 

161.0 


2010 May 1-3, 20-24, Jun 1, 5-9, 11 14, 27 

6C 

40.5 

ELAIS-S1 C1241 J>,c ’ d 

2004 Jan 9-11 

6A 

24.6 


2004 Jan 30, Feb 1 

6B 

18.6 


2004 Dec 19, 27, 2005 Jan 1-3 

1.5D 

40.2 


2005 Jan 9-11, 20-22 

750B 

50.0 


2005 Mar 25, Apr 8, 11 

6A 

27.2 


2005 Apr 24, 26, 30, May 1 

750A 

34.3 


2005 June 8, 9 

EW367 

18.3 


2005 June 19, 24 

6B 

18.5 

C1967 d 

2009 Dec 19-31, 2010 Jan 1-3, Apr 1-6, 12, 14-18 

6A 

68.6 


2009 Aug 13, 19, 21-27 

6D 

87.9 


2009 Nov 10-18, 20-29 

6B 

9.8 


2010 Apr 19 

6C 

10.6 


a Data presented in i Norris i ^tjih (2006). 
b Data presented in ^iddelberg^t^L (2008) . 
c Data presented in JHales^etjiL (2014a). 
d Data presented in this work. 









ATLAS Third Data Release 


19 


Table A2. Coordinates of mosaic pointings in the two ATLAS fields. 


Pointing 

R.A. (J2000) 

Decl. (J2000) 

CDFS 



1 

03:28:47.33 

-28:38:37.98 

la 

03:27:18.36 

-28:38:31.14 

2 

03:28:03.89 

-28:21:46.74 

3 

03:28:48.48 

-28:05:05.58 

3a 

03:27:18.36 

-28:05:05.58 

4 

03:28:05.26 

-27:48:14.34 

5 

03:28:49.61 

-27:31:32.82 

5a 

03:27:18.36 

-27:31:32.82 

10 

03:30:16.97 

-27:31:40.02 

11 

03:29:32.83 

-27:48:22.98 

12 

03:30:16.30 

-28:05:12.42 

13 

03:29:31.92 

-28:21:55.74 

14 

03:30:15.60 

-28:38:44.82 

15 

03:31:43.87 

-28:38:48.42 

16 

03:30:59.95 

-28:22:00.78 

27 

03:32:27.99 

-28:22:02.58 

28 

03:33:12.12 

-28:38:48.42 

29 

03:34:40.39 

-28:38:44.82 

30 

03:33:56.02 

-28:22:00.78 

31 

03:34:39.70 

-28:05:12.42 

33 

03:34:39.03 

-27:31:40.02 

41 

03:32:28.00 

-27:48:30.00 

42 

03:31:20.17 

-27:48:30.00 

43 

03:31:54.08 

-28:01:29.44 

44 

03:33:01.92 

-28:01:29.44 

45 

03:33:35.83 

-27:48:30.00 

46 

03:33:01.92 

-27:35:30.56 

47 

03:31:54.08 

-27:35:30.56 

ELAIS-S1 



1 

00:32:03.55 

-43:44:51.24 

2 

00:31:10.95 

-43:27:59.64 

3 

00:32:05.04 

-43:11:18.84 

4 

00:33:51.29 

-43:11:24.96 

5 

00:32:57.67 

-43:28:09.00 

6 

00:33:50.79 

-43:44:57.36 

7 

00:35:38.02 

-43:44:57.36 

8 

00:34:44.40 

-43:28:11.88 

9 

00:35:37.51 

-43:11:24.96 

10 

00:37:23.76 

-43:11:18.84 

11 

00:36:31.13 

-43:28:09.00 

12 

00:37:25.25 

-43:44:51.24 

13 

00:36:31.13 

-44:01:42.84 

14 

00:37:25.25 

-44:18:34.44 

15 

00:35:38.02 

-44:18:34.44 

16 

00:34:44.40 

-44:01:42.84 

17 

00:32:57.67 

-44:01:42.84 

18 

00:33:50.79 

-44:18:34.44 

19 

00:32:03.55 

-44:18:34.44 

20 

00:31:10.95 

-44:01:42.84 






